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I. INTRODUCTION 



Since the original work of Shirkov and Solovtsov [1] appeared in 1997, the analytic ap- 
proach to QCD perturbation theory has evolved and progressed considerably. At the heart 
of this approach is the spectral density which provides the means to define an analytic run- 
ning coupling in the Euclidean space via a dispersion relation in accordance with causality 
and renormalization-group (RG) invariance. Using the same spectral density one can also 
define the running coupling in Minkowski space by employing the dispersion relation for the 
Adler function [2-5]. In parallel, this approach has been extended beyond the one- loop level 
[5, 6] and analytic and numerical tools for its application have been developed [7-12]. These 
efforts culminated in a systematic calculational framework, termed Analytic Perturbation 
Theory (APT), recently reviewed in [13]. 

The simple analytization of the running coupling and its integer powers has been gener- 
alized to the level of hadronic amplitudes [14, 15] as a whole and new techniques have been 
developed to deal with more than one (large) momentum scale [16-18] (for a brief exposition, 
see [19]). This encompassing version of analytization includes all terms that may contribute 
to the spectral density, i.e., affect the discontinuity across the cut along the negative real 
axis — oo < Q 2 < 0. It turns out that logarithms of the second large scale, which can be 
the factorization scale in higher-order perturbative calculations or the evolution scale, cor- 
respond to non-integer indices of the analytic couplings, giving rise — in Euclidean space — to 
Fractional Analytic Perturbation Theory (FAPT) [20, 21]. This concept was successfully 
extended to the timelike region and a unified description in the whole complex momentum 
space was achieved [22] (see also [23, 24]). The issue of crossing heavy-flavor thresholds, 
naturally entering such calculations, has been considered within APT [13, 18, 25, 26] and 
more recently also within FAPT [27]. 

Another important topic, which is at the core of the present investigation, deals with the 
perturbative summation in (F)APT. As it has been demonstrated in [28-30], it is possible to 
determine the total sum of the perturbative expansion at the level of the one-loop running 
of the coupling. This result will be extended here to include in the sum at least the two-loop 
running of the coupling. Even more, making a natural assumption concerning the asymptotic 
behavior of the perturbative coefficients, proposed long ago by Lipatov [31], we will show 
that it is possible to estimate the sum of the series to all orders of the expansion. This 
important feature of the (F)APT non-power series gives us the possibility to estimate the 
uncertainties of the perturbative results with higher accuracy relative to the conventional 
QCD power-series expansion. In the present investigation we will extend and systematize 
these issues towards a complete calculational scheme considering also some applications- 
including updated predictions for the decay width of a Higgs boson to a bottom-antibottom 
pair, relevant for the Higgs-boson search at the Tevatron and the LHC. 

The paper is organized as follows. In Sec. II we compare the results obtained in different 
schemes— FOPT, CIPT, and FAPT— and prove that CIPT and FAPT provide identical 
results for -R e + e -_>hadrons an d Rs- 111 addition, we recall in this section the pivotal ingredients 
of the analytic perturbative approach and discuss briefly its extension to the case of a variable 
flavor number (the so-called "global case" [25]). Moreover, in Sec. Ill we describe how 
typical series expansions within FAPT can be resummed at the level of the one-loop coupling 
running having recourse to a generating function — referring the reader for the two-loop case 
to Appendix D. Section IV is devoted to the Higgs-boson decay into a bb pair, a process which 
contains and exhibits the conceptual and technical details of our analytic framework. In the 
same section, we calculate the associated decay width, which involves the Adler function of 
a scalar correlator, and compare it with the results of the standard perturbative scheme. 
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Then, in Sec. V, we turn to another application of the proposed resummation technique and 
consider the Adler function of a vector correlator, pertaining to -R e + e -_j.hadrons- Finally, Sec. 
VI contains our concluding remarks, while some important technical derivations are given 
in five Appendices. 



II. FOPT, CIPT AND FAPT 

A. Two-point correlator of scalar/vector currents 

As mentioned in the Introduction, the initial motivation to invent new QCD couplings 
was the desire to interrelate the Adler D-function, 

D{Q\ /i 2 ) ^ D{Q 2 ) = d + J2d m aT(Q 2 ) , (1) 

m=l 

a(e + e~ — > hadrons) 

calculable m the Euclidean domain, and the quantity R e + e - = — , 

a(e + e — > ) 

R ( s , ^ R ( s ) = ro + J2 r m aT(s) , (2) 



m=l 



which is measured in the Minkowski region. Both quantities are considered in standard QCD 
perturbation theory, demanding that the couplings entering them satisfy the renormalization- 
group equation. To facilitate direct comparison of our results further below with the higher- 
order calculations in standard perturbation theory, in particular Refs. [32-34], we use here 
the variable a s = a s /n. Thus, the beta-function coefficients of this coupling are defined by 



A(fls) = rfln(/i 2 ) = ~ Us ^° a s + Pi a s + ■■■) = - — 



2 

+ . . . 



(3) 



where /3 n = b n /A n+1 and the coefficients b n are specified in Appendix D. 

The functions D and R can be related to each other via the following dispersion relations 
without any reference to perturbation theory 

fl « = k ^ 27, £ 1*'^ d ° • D ^ ~ d\R\ = Q 2 [ j^d, , (4) 

where for the first term in Eq. (4), the integration contour I\ around the cut (solid red 
line) is shown in Fig. 1. However, employing a perturbative expansion on the LHS of Eqs. 
(1) and (2), one obtains, in fact, a relation between the powers of ln(s//i 2 ) and \n(Q 2 /fi 2 ) 
in the coefficients r m (s,/i 2 ) and d n (Q 2 ,fi 2 ), while the powers of the couplings a s (/x 2 ) reveal 
themselves as parameters. 



B. Fixed-order perturbation theory 



For the fixed order-n perturbation theory (abbreviated as FOPT), one can start from Eq. 
(4) and use (2) on its LHS and (1) on its RHS to obtain 



R, 



FOPT, 



s) = R n (s,fi 2 = s) 



1 

2vn 



-s+ie 



D n (a, /i 2 



■ da 



r + ^2r m (a s (s)Y 



(5) 



fl 2 =S 



m=l 
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To further utilize the FOPT approach, it is useful to consider the relation between the 
coefficients r m and df. in more detail. The goal is to express the coefficients r m in (5) in terms 
of the calculable coefficients d n in Eq. (1), i.e., to write r m = T mk dk, where summation over 
k = 1, . . . , m is implied. The matrix T mk is triangular with unity elements on its diagonal — 
see Table I. In the horizontal direction, i.e., along the rows of this Table, we include all 
coefficients di up to the coefficient d 5 , the latter not calculated yet, but due to be estimated 
within our approach later in connection with specific applications — Sees. IV and V). 

The elements, proportional to which originate from the one-loop evolution procedure, 
have the following general form 



<^ (n-l + 2m)! (a s 7r^ 
n 4-^- (n-1)! (2m)! 1 ' 2m + 1 1 



m=0 



and can be obtained for any fixed order n of the expansion by the procedure described in 
Appendix A. The other /^-coefficients — related to higher loops — have been color-printed 
below using the same color assignments as in Table I. 

To get acquaintance with the use of this Table, we write out explicitly the relation between 



TABLE I: This table exemplifies the structure of the first few coefficients r m = T dk (summation 
over k = l,...,m implied) of the conventional expansion of the -R-ratio, where T mk are the table 
entries. Each coefficient r m contains a number of dk (k < m) terms in its expansion that are shown 
in the corresponding row. The various contributions are marked by different colors according to 
their loop order: one-loop — black; two-loop — red; three-loop — blue. 
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"This expression for r$ was also obtained in Refs. [35, 36]. 
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(7) 



The role of the "kinematical ir 2 " -terms becomes more pronounced for higher orders m. For 
instance, as regards the term r 3 , the underlined ^-contributions are comparable in size with 
the original G?3-contribution, while for r 4 these contributions even exceed in magnitude the 
value of the coefficient d\ [33, 34], as one observes from the expressions 



r 3 = 18.2 - 24.9 + (-4.22 + 3.02)71/ + (-0.086 + 0.091 )n : 
r 4 = 135.8 - 292.4 + (-34.4 + 53.2)71/ + (1.88 - 2.67W 



/' 



/ 



(8) 

-0.01 + 0.032 )n 3 f . (9) 



Therefore, we may be tempted to take into account the "kinematical" 7r 2 -terms to all orders 
of the expansion, i.e., to sum, for a fixed element d n , over the index m along a single column 
in Table I by taking into account the corresponding power of the coupling a™. 
Let us now look at alternative QCD perturbative expansions. 



C. Contour-improved perturbation theory 

Another way to determine R was suggested in [37, 38], where the integration contour in 
Eq. (4) was changed from the cut Ti to the contour I^, i.e., to the circle in the complex Q 2 - 
plane with radius s centered at 0, see Fig. 1. Applying the FOPT expansion, the integration 
of the terms d n (a/fi 2 ) a"(/J 2 ) in Eq. (1), along the contour T 2 , is completely equivalent to 
the integration along the cut. Using instead the so-called Contour-Improved Perturbation 
Theory (CIPT), for large enough values of s, one can integrate the terms d n a^(a), entering 
its RHS, along the contour T 2 to get 

A* = <{a)-. (10) 

2-ni J F2 a 

This means that one may absorb all logarithms inside the running coupling just by adjusting 
in Eq. (4) the magnitude of the scale fi 2 = a — without performing any expansion. CIPT 
is currently considered to be the most preferable technique to account for the running of 
perturbative quantities in a number of processes, including the r-decay [35, 39], the R e + e -- 
ratio [33, 34], and the width of the Higgs-boson decay H— > bb [40]. The integration of the 
running couplings along the contour T 2 is performed numerically [35], but for the one-loop 
running the result for i? CIPT can be obtained also analytically and is found to coincide with 
expression (6), as expected (see Appendix A). 



The four- loop-order contribution has also been computed in Ref. [36]. The two results coincide. 
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Re z = Q 2 





FIG. 1: Integration contour T\ entering the calculation of R(s) in Eq. (4). The position of the 
Landau pole is marked by the cross symbol X . The other integration contours shown are explained 
in the text. 

D. Fractional analytic perturbation theory 

Inspection of Table I suggests to consider the sum of the elements of each of its columns 
as defining a new coupling. For instance, the first column, associated with d\, gives rise to 
the coupling 2ti, while the second column pertains to the coupling 2t 2 , and so on. In this 
way, one becomes able to introduce a non-power-series expansion in Minkowski space. This 
is, actually, just another way to define the Analytic Perturbation Theory [1, 4, 13], and its 
extension — Fractional Analytic Perturbation Theory [20, 22] (see also [14, 21] and [24, 27] 
for reviews). The basis of (F)APT is provided by the following linear operations which define 
analytic images of the normalized coupling a = a s 6 /4 and its powers in the Euclidean and 
the Minkowski space, respectively: 2 



where p u (cr) = — Im [a u (— a)] is the spectral density. The set of the couplings {A u , %l u } 

7T 

satisfy the dispersion relations Eq. (4) and fulfill the constraint RD = 1. Applying the 
operation Am on the RHS of the last expression in Eq. (1), one obtains R in APT: 



2 To streamline our notation, we use for all quantities in the Euclidean space a calligraphic symbol, whereas 
their Minkowski-space counterparts are denoted by Gothic symbols. Moreover, note that Latin indices 
represent integer numbers n, whereas Greek indices v mark real numbers. The subscripts E and M serve 
to specify the space we are working in: Euclidean and Minkowski, respectively. 




1 





m=l 
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The reader should note here that the expansion coefficients d n differ from those in Eq. (1) 
and are defined as follows 



d i, 



(13) 



The prefactors (4/6 ) n above — and analogously in Eq. (14) below — serve to connect this 
analysis to our standard definitions of the analytic couplings A v and 21^ in Refs. [20, 22]. 
Hence, according to Eqs. (11)— (12), one may associate the evaluation of R within (F)APT 
with the integration contour T 3 in Fig. 1. 

Consider now how the (F)APT result R APT in Eq. (12) correlates with the elements of 
Table I. One appreciates that every term d m 2l m of the series appears as an infinite sum of 
the elements along the corresponding matrix column in this table. The series is convergent 
and includes all so-called "kinematical terms", as we discussed at length in [22]. One can 
verify this by considering the expansion of the first few elements d m 2l m in Eq. (14) and then 
compare the results with the content of the corresponding column in Table I. A further 
advantage of Eq. (14) is that one can use it to derive expression (6) in a more direct way, 
namely, as an expansion of the one-loop Minkowski analytic couplings 2l n in powers of 
the variable a s — terms printed in black color in Eq. (14). The analogous terms for higher 
loops can be found in Appendix C of Ref. [22]. Thus, we obtain (with the same coloring 
assignments as in Sec. II B) 



(3) 



1 - 



(a s 7r/3 



+ 



(a s 7r/3 ) 



7T 2 , / 5 „3 

— at (/3 /3 1 - + /3j- a 



+ 



7T 2 „ 77 2 
— PlP!— a 2 + 
5 4 



7T 



-—all f3 f3 2 3 + fcfo -a s + 



7T 



a 6 s /3 (3 3 7 + ... 



(14a) 



- a 



(2) 



[1 - (a s vr/3 ) 2 + (a s vr/3 ) 4 . . .] 



7T 



(/3 /3i7 + /3 2 4a s + ...) 



7T 



- a 



(2) 



a 6 s f3 f3 2 8 + ... 

(a s 7r/3 ) 2 „ 



(14b) 



7T 2 _ 27 

-— atPofc— + 



(14c) 



21 



(2) 



1 - 



{a 2 s n(3 ) 



■10 + . 



+ 



(14d) 



The main conclusion from the above exposition is that, depending on the particular scheme 
of the perturbative expansion used, the elements of Table I may be summed in different ways 
(see Table II). Specifically, 
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(i) FOPT gives the sum of a finite number of terms along some row to create r m , and 
then — following Eq. (5) — it sums the results up to n to yield R n , i.e., R n = d + 

n 

E a™T mk d k . 

m,k=l 

(ii) (F)APT takes into account each infinite column as a whole in the form of the expansion 
J m 2l m , thus including this way all "kinematical terms" by construction, and then sums 
a number of 2l m terms into -R^ PT in the form given by expression (12). 

It is evident that for any fixed order n, the results for R n and Rn FT cannot coincide. Nev- 
ertheless, it was shown in [22] that calculating the decay width of a Higgs boson into a bb 
pair using CIPT, leads to the same result one would obtain for R using FAPT with one-loop 
running of the coupling. This coincidence turns out to be not accidental but to hold at any 
loop order of the perturbative expansion by virtue of the dispersion relation given in Eq. 
(11), as we will prove next. 

TABLE II: QCD observable R calculated with different perturbative expansions: (F)APT — 
(Fractional) Analytic Perturbation Theory; CIPT — Contour-Improved Perturbation Theory; 
FOPT — Fixed-Order Perturbation Theory. The associated contours T\, T2, T3, are displayed in 
Fig. 1. 



Perturbative scheme 


Contours 


R expansion 


FOPT 




E„ r n< 


CIPT 


r 2 


d n A n 


(F)APT a 


F 3 





a The coefficients d n and d n in CIPT and (F)APT differ by trivial factors, see Eq. (13), due to the different 
normalization of the couplings A n and 2(„. 



E. Fixed flavor number vs. global FAPT 

We commence our analysis within FAPT by recalling the salient features of the analytic 
approach to QCD perturbation theory, expanding our remarks given in Section II D. In 
order to have a direct connection to our previous papers on the subject [20-22], and to 
simplify the main formulae in those sections where we consider fixed-order (F)APT with a 
constant value of active flavors Nf, we use here the normalized coupling of perturbative QCD 
(pQCD) [28] a = (3 f a s with the useful abbreviation (3 f = b (N f )/(4ir) = (11-2 N f /3)/(Air), 
where bo(Nf) denotes the first coefficient of the QCD (3 function. Analytic images of the 
normalized coupling and its powers are constructed by means of the linear operations A E 

and A M according to (11) using the spectral density p u (a) = — Im [a v {— a)\. 

To be in line with the above definitions, we also introduce analogous expressions for the 
fixed- Nf quantities with standard normalization, i.e., 

A(Q 2 ;^) = ^p, Ms;N f ) = ^, (is) 

which correspond to the analytic couplings A v and 2l„ in the Shirkov-Solovtsov terminology 
[13]. These couplings have dispersive representations of the type (11) with spectral densities 
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p u (cr;Nf) = p v (cr)/f3f. For the sake of simplicity we will omit to display Nf (and other 
evident arguments) explicitly In the present analysis we express all variables in terms of 
L = ln(Q 2 /A 2 ) (Euclidean space) or L = ln(s/A 2 ) (Minkowski space), using the notation 
defined above, but with an argument placed in square brackets: a u [L], A U [L], and 2l„[L]. 
Then, in the one-loop approximation (labeled by the superscript (1)), we have 



— , Sli 1 ' 1 [L] = — arccos 

'i L — 1 7T 



L 



y/L 2 + 7T 2 



P?[L;N f 



L 2 + 7T 2 



7T Z 



(16a) 
(16b) 



On the other hand, when we will discuss the global version of (F)APT, where Q 2 (or s) varies 
in the whole ("global") domain [0, oo), and Nf effectively becomes dependent on Q 2 (or s), 
we will use in Eq. (11) the global version of the spectral density that takes into account 
threshold effects and is, therefore, Aff -dependent. 

In order to make the effect of crossing a heavy-quark threshold more plausible, consider 
a single threshold at a 



(corresponding to the transition Nf = 3 =>■ Nf = 4) and write 
the spectral density in a form which expresses it in terms of the fixed-flavor spectral densities 
with 3 and 4 flavors, p n [L; 3] and p n [L + A4; 4], respectively. The result is 



pt° h [L*\ 



9 (L a < L A ) p v [L a ; 3] + 9 (L 4 < L a ) p v [L a + A 4 ; 4] 



(17) 



with L a = In (cr/A 2 ), L4 = In (ml/ A 2 ), and A4 = In (A|/A 2 ). Note, however, that in standard 
QCD perturbation theory the coupling at the threshold becomes discontinuous. To enable 
a smooth matching at this point, one has to readjust the values of Aq CD on each side of the 



threshold in correspondence with their associated flavor numbers (Aj = A 



QCD 



In Fig. 2 



we show the global spectral density pf lob [^] in comparison with the spectral density pi[L;3] 
which corresponds to a fixed flavor number Nf = 3. We see that, because of the different 
values of A 3 and A 4 , and also of ^3 and ^4, these quantities differ from each other when L 
crosses L 4 . 

Substituting the obtained spectral density pf oh (a) [cf. Eq. (17)] into Eq. (11), we get 
continuous expressions for the analytic couplings in both domains of the complex Q 2 space. 
In the Minkowski region, the global analytic coupling reads 

2lf h [L] = 6 (L<L 4 ) 3] - % [L 4 ; 3] + %{L A + A 4 ; 4]} + 6 (L 4 <L) %[L + A 4 ; 4] , (18) 



0.1 



. 04 



. 02 




2 2 4 



FIG. 2: Spectral densities pf lob [-L] (solid line) and pi[L;3] (dashed line). The discontinuity at 
L = L 4 is artificially enhanced by a factor of 1.5 to make it more visible. 
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whereas its Euclidean counterpart assumes the form (referring for more details to [27, 29]) 

U _ 

AT[L] =A U [L + A 4 ; 4] + f PAL^] -p [L +X 4 ;A} ^ 

J J. I a & 

oo 

To demonstrate the magnitude of the threshold corrections, we show in Fig. 3 the values of the 
normalized deviations AA U [L] = A s J ob [L] -A u [L+\ 4 ; 4] and A2f„[L] = 2t;f b [L] -&„[!,+ A 4 ; 4] 
in the Euclidean and the Minkowski domain, respectively. We see that in both domains 




-10 -5 5 10 -10 -5 5 10 



FIG. 3: Left: Deviation of the global coupling relative to the fixed-iVy coupling in FAPT: 
AAl [L]/Af oh [L]. Right: The same for A2li[L]/2tf b [L]. 

these deviations vary from —20% for large values of — L w 10, going through zero in the 
vicinity of L « —5, and then increase up to the value +20% for L w 0, tending, finally, to 
as L — > oo. This means that these deviations reach the level of 10% in the region of several 
tens of GeV 2 . 



F. Relation between CIPT and (F)APT 

To establish the equivalence between FAPT and CIPT, we consider a more general expan- 
sion than Eq. (1) which contains the coupling with a non-integer power that can be related 
to an anomalous dimension. Such a quantity reads 

D v = d a: + J2d n a u s +n , (20) 

n=l 

where v is not an integer number. 

Symbolically, we have the following equivalence 

FAPT ^ J = CIPT (/ ) ' ( 21 ) 

referring for the integration contours to Fig. 1. To establish the equivalence between FAPT 
and CIPT in the above relation, we employ the chain of the equalities 

FAPT . r P -^l da s rim [„*(-„)] j, = ± f da 

J s a J s vr a 2mJ r3 K J a > 1 } 



10 



where in the last step the integral has to be evaluated along the contour T3. On the other 
hand, the CIPT part of Eq. (21) can be identically rewritten as 



CIPT = -L [ a»^. 
2m J T2 y ' a 



(23) 



Therefore, we have to prove that 




v da 



+ 




V 



da 



0. 



(24) 



Jr 3 



a 



Jv 2 



a 



To this end, we close the contour T 3 along the large circle T 4 with radius R that tends to 
00, and take into account the closed composed contour r 234 = T 2 U T 3 U T 4 (see Fig. 1). The 
integral /(r 234 ), 



along the closed contour r 234 is equal to the sum of the residues inside the enclosed region. 
Provided the radius s of the contour r 2 is large enough, no poles from a(a) owing to per- 
turbation theory are inside the contour r 2 3 4 . As a result, 7(^34) = 0. Moreover, (a(<r)) I/ 
(u > 0) on the contour T 4 decreases with growing R, and therefore this contribution vanishes 
as R — ?■ 00. Consequently, we finally obtain 



The above equivalence not withstanding, there is a crucial advantage of the FAPT ap- 
proach with respect to CIPT. In fact, as long as one is only interested in a numerical 
estimate, CIPT provides acceptable results for several typical processes. However, if one 
pretends to employ analytic expressions, and thereby to control each step of the calculation, 
CIPT is not sufficient. In that case, one needs another perturbative scheme that is able to 
yield explicit expressions for the couplings along each column in Table I. Such a scheme is 
naturally provided by (F)APT. Moreover, we shall show below that this scheme can even 
admit the resummation over columns [29]. This is an important feature, given that the 
conventional perturbative series in Eq. (1) cannot amount to a unique, i.e., resummation- 
method-independent, result owing to the asymptotic nature of the power series. 

III. RESUMMATION IN (F)APT IN THE ONE-LOOP ORDER 

In this section, we consider different sorts of perturbation-series expansions of typical 
physical quantities, like the Adler function, D[L\. Our goal is to perform the summation of 
such expansions under the imposition of a couple of basic constraints. As we will show in 
a moment, these constraints are (i) a recurrence relation for higher-order couplings and (ii) 
the nonpower character of the series. These considerations will be based on an appropriate 
generating function for the expansion coefficients. 

In what follows we discuss in detail the one-loop running case. However, also the techni- 
cally more complicated two-loop running case is worked out and the corresponding expres- 
sions are provided in Appendix D. 




(25) 




(26) 
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A. Generating function for the series expansion 



At the one- loop level [22], we have 

D[L] 




71=1 




do + d^dn { ML] \ (27) 



with d n = d n /di [cf. Eq. (13)], where L = ln(Q 2 /A 2 ) applies to the Euclidean quantities 
(D[L], D APT [L], a n [L], and „4 n [L]), and L = ln(s/A 2 ) pertains to their Minkowski versions 
R APT [L] and 2l n [L]. 

The resummation of this series is on the focus of the present work. For this purpose, it is 
useful to introduce a generating function P(t) for the series expansion and write 



d, 



poo poo 

= / P(t) t^dt with / P(t) dt = \. (28) 
Jo Jo 



We will show in the next step how to use this generating function in order to resum nonpower 
series expansions (both in the Euclidean and the Minkowski region). But first recall that 
the standard coupling a n of perturbative QCD, as well as the couplings „4 n ,2l n — together 
with the spectral density p n — satisfy a one-loop renormalization-group equation that can be 
recast in the form of a recurrence relation [41]: 



a n+l [L] \ _ 1 / d \ n / a^L] 



F n+1 [L\ J r(n + l) \ dLj { Ti[L] 



(29) 



Here ^[L] denotes one of the analytic quantities A[L], 2l[L], p[L]. Substituting Eqs. (28) 
and (29) into the perturbative-series expansion, i.e., into Eq. (27), one obtains [28] in analogy 
to the previous equation, 

where S = D APT or R APT — depending on the choice of the analytic couplings J 7 = A or 
T = 21, respectively. In the above equation, and in the considerations to follow, we use the 
abbreviated notation 

POO 

{{T[L - t])) p = T\L — t]P(t) dt . (31) 



As long as we have not proved that summation and integration can be interchanged, this 
representation has only a formal meaning. Note, however, that the integration over the 
Taylor-series expansion of the term a[L — t] in the integrand reproduces the initial series for 
any partial sum. The integrand in the standard case of the QCD running coupling (first 
line in Eq. (30)) faces a pole singularity (termed the infrared renormalon singularity) and is, 
therefore, ill-defined. In contrast, the integral in the second entry in Eq. (30) has a rigorous 
meaning by virtue of the finiteness of J"i[£], being one of the quantities {«4i[L], 2li[L] < 
1) Pi[L] < I/ 71 " 2 }- Therefore, the expression on the RHS of Eq. (30), together with (31), can 
be called the sum of the corresponding series in the sense of Euler. Since any coefficient d n 
is the moment of P(t) [cf. Eq. (28)], this function should fall off faster than any power — e.g., 
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like an exponential or faster. 3 Therefore, all APT expressions on the RHS of Eq. (30) (second 
line) exist and are proportional to T\L — i(L)], where, for each L, t(L) is the average value 
of t associated with this quantity. 

Provided the generating function is known, one can compute the integral (31) explicitly 
and obtain this way all-order estimates for the expanded quantity for any series expansion, 
provided the coupling parameters fulfill the following conditions: 

(i) satisfy the one-loop renormalization-group equation [cf. Eq. (29)], 

(ii) are real, and 

(iii) have only integrable singularities. 

If the first coefficient do of the expansion (27) is not accompanied by unity (a ), but by some 
other fractional power v of the coupling (a u ), then, as it has been shown in [27, 29], the 
resummation method (30) has to be modified to read 

S V [L] = d T V [L] + d x {(F 1+V [L - t])) Pv , (32) 

where now the generating function P v depends also on v. This quantity can be deduced 
from P as follows: 

Pu{t)= f'pf-L-)^)-^-. (33) 
Jo \l-xj l-x 

Here $ v (x) = vx v ~ x , so that lim^o &v — > $( x ), an d therefore lim^o Pvif) = P(t)- The last 
step completes the generalization of the original APT resummation procedure of [28] to the 
case of FAPT. 



B. Modeling the expansion coefficients and their generating function 

For most relevant QCD processes, only the first few coefficients d n are known, while 
the computation of higher-order coefficients is technically a highly complicated task. This is 
despite the impressing development of sophisticated algorithms during the last few years [32- 
34, 43, 44]. In view of this, it is extremely useful to have alternative methods for calculating 
the higher-order coefficients, let alone to resum the whole series — even if the result represents 
a sort of approximation — provided the quality of the applied method is high and the inherent 
uncertainties entailed can be kept under control. On the other hand, the asymptotic form 
of the coefficients d n can be predicted from 

d n ~ T(n + 1) tfc 11 [1 + 0{l/n)] -> P(t) ~ {t/ C y i+1 e- t/c (34) 

a form that is inspired by Lipatov's asymptotic expression in Ref. [31] (see also [45], and [46] 
for a review), and where 7 < 1 and c are numerical coefficients. One anticipates that the 
large-order behavior of the expansion coefficients translates into the asymptotic form of the 
generating function P(t). 

To proceed, we have to construct first a model for the expansion coefficients d n = d n /di, 
having recourse to the information about the first few fixed-order coefficients. The second 



3 Let us mention in this context that an expression similar to (30) was also obtained in [30] . The authors of 
this paper have used the "large /?o" approximation to create a specific model for the generating function 
in conjunction with Neubert's approach [42]. 
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step is to interpolate between the obtained result and the expression for the tail obtained for 
d(n) from (34) at asymptotically large orders n. For a fixed-sign series we adopt the model 
(note the normalization d™ del = 1) 

-del Wd + A/-^^-!)! 

dn = a^Ta 2 ' (35) 

where the parameters 7^, and c are determined from the values of the first few known 
coefficients d n . We mention incidentally that the case of an alternating-sign series has been 
considered in [28] and will not be addressed here. As regards the most general case of the 
series, the behavior of <i™ odel with respect to n turns out to be the same as what is obtained in 
the renormalon approach (that is usually supplied with the Naive Non-Abelization (NNA) 
[47]), see, e.g., [48] — except perhaps for the specific values of the parameters which are 
different. Modeling the expansion coefficients according to Eq. (35), leads to the following 
generating function P(t): 

pmodel^ _ [Ax(t/ C ) 1+71 + A 2 (t/ C ) 72 ] 6~ t/c . (36) 

Using the first already known coefficients, we can calibrate our model in order to extract 
information about still higher and unknown orders, as well as to gain information on the 
resumed behavior of the whole series. It should be emphasized that for large values of 
the argument t, our model for P(t) can become rather crude. This, however, does not 
significantly change the final result of the summation due to the convergence of the integral 
in Eq. (31). For this reason, it is not necessary to know the asymptotic behavior of p model (£) 
very accurately. All said, let us now consider concrete examples to understand the modus 
operandi and the benefits of this technology. 

IV. MASTER EXAMPLE: HIGGS-BOSON DECAY INTO A lb PAIR 

Our goal in this section is the calculation of the width of the Higgs-boson decay into a bb 
pair, i.e., 



rw b (M„) = r^K 2 ) , (37) 

3m b 

where r^m 2 ) = 3G F M H m 2 b /4:V2ir, m b and M H are the pole mass of the 6-quark and the 
mass of the Higgs boson, respectively, and R^Mjj) = m 2 (M^) i? g (M^). Our interest in 
this quantity derives from the fact that this process contains all main ingredients for such 
high-order calculations, discussed above, with known coefficients up to the order 0(aj) [32]. 
The quantity Rs(s) is the discontinuity i?s(s) = Imll(-s — ie)/(2n s) of the imaginary part 
of the correlator II of the two scalar currents Jjf = ^b^b for bottom quarks with an on-shell 
mass mb coupled to the scalar Higgs boson with mass Mh, and it is given by 

H(Q 2 ) = (4tt) 2 * J dxe^(0\ T[ J fe S 0rK S (O)] |0> , 

where Q 2 = — q 2 . Direct multi-loop calculations are usually performed in the far Euclidean 
(spacelike) region for the corresponding Adler function D$ [32, 40, 49, 50], where QCD 
perturbation theory works reliably. Hence, we write 



q 2 d r u{Q 2 

6 dQ 2 V Q 2 



D S {Q* 



l + J2d n a:(Q 2 ) 



n>l 



(38) 
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A. Generating function for the scalar correlator and estimates for higher-order 
coefficients 



The Adler function, related to the scalar correlator (see [22, 32]) and pertaining to the 
Higgs-boson decay, reads 

oo 

D S [L] = 1 + Y^dn (a s [L]) n . (39) 

n=l 

Note that here the definition of the coefficients d n is given by d n = d n jd\. The n-dependence 
of the coefficients d n , in accordance with expression (35), can be simulated by the following 
two-parameter model 

% = <T- 1 ^r(n), (40a) 

a form which ensures that df (the superscript H denoting "Higgs") is automatically equal 
to unity. Concurrently, these coefficients can be reproduced by the following generating 
function: 

« 40b » 

To test these formulae, we use the two known coefficients d 2 and d 3 for this process and 
plug them into model (40a) aiming to predict the next two coefficients df and df. The 
result for the next coefficient df is shown in the second row of Table III: df = 662. We see 
that it is fairly close to the value 620 calculated recently by Chetyrkin et al. in [32]. This 
procedure can be geared up to predict the next coefficient d$. Indeed, taking into account the 
coefficient d^ and slightly readjusting the parameters c and 5 of our model (40b) from their 
previous values {c = 2.5, 5 = —0.48} to {c = 2.4, 5 = —0.52}, we find (third row in Table 
III) d 5 PS 7826. This result agrees very well with the value 7782, derived in [36] by appealing 
to the Principle of Minimal Sensitivity (PMS). Moreover, it has a reasonable agreement with 
our previous models (rows 3 and 4 in Table III) and also with the improved NNA (INNA) 
prediction shown in row 7, while the original NNA prediction (value given in row 6) fails. We 
may conclude that our model calculation of the coefficient d 5 is congruent with the prediction 
following from the application of the PMS method. Both these approaches provide results 
for most coefficients d n approximately twice larger than those found by applying the INNA 
approximation. 



B. Higgs-boson decay width in FAPT 

Within the one- loop approximation of FAPT, the quantity R S (M^) has the following 
non-power series expansion: 4 

i?r T [L] = 3 1 gl ° b M + d i E % W } • ( 41 ) 



The appearance in the denominators of the factors 7r™ in tandem with the coefficients d n is a consequence 
of the particular d n normalization — see Eq. (39). 
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TABLE III: Perturbation Theory (PT) coefficients d n for the Z?s series obtained from different 
calculations. 





PT coefficients 


di d 2 d 3 d 4 d 5 


1 


pQCD results from [32] 


1 7.42 62.3 620 




Models and estimates 




2 


Model (40) with c = 2.5, S = -0.48 


1 7.42 62.3 662 8615 


3 


Model (40) with c = 2.4, 5 = -0.52 


1 7.50 61.1 625 7826 


4 


"PMS" predictions from [32, 36, 51] a 


- - 64.8 547 7782 


5 


"NNA" predictions from [50, 52] 


1 3.87 21.7 122 1200 


6 


"INN A" predictions from App. C b 


1 3.87 36 251 4930 



°The value of the coefficient d$ has been computed here with the method of [36]. 
b We employ the upper limit of the INNA predictions — consult Appendix C. 



Here is the renormalization-group invariant mass satisfying the one-loop m 2 (/i 2 ) evolu- 
tion equation 

mfaapffl (42) 

with vq = 27o/6o(5) = 1.04, where 70 is the quark- mass anomalous dimension, and the 
renormalization-group invariant quantity TO(i) is denned through the effective RG mass 
m b (ml). It is worth recalling that Eq. (41) in the one-loop running of the coupling coincides 
with the result obtained in [50] for R s using CIPT. 

The key question now is to which extent FAPT is able to reproduce with a sufficient 
quality the whole sum of the series expansion of i?g APT [L] in the range L e [11.7, 13.6]. 
This will be checked by taking recourse to the model (40) in conjunction with Eq. (B7). 
Note that this L range corresponds to the Higgs-mass values M H e [80, 180] GeV with 

Aq^ 3 = 201 MeV and 2li° b (mi) = 0.1226. In this mass range, we have L 5 < L < L 6 , so 
that Eq. (B7) transforms into 



(43) 



with P UQ (t) (defined via Eq. (33)) and where we have evaluated Eq. (40) with the parameter 
values c = 2.4, 5 = -0.52, and v = 1.04. 

The all-order expression (43) above allows us to determine the accuracy of the truncation 
procedure of the FAPT expression 



N 7 



rS(i; n] = r^(m 2 1} ) |s^[l] + 4^-^ a*£jL] \ (44) 

at order and estimate the relative errors 



n=l 



A n[L] = 1 - ^i 1 '^ ■ (45) 

H^bb 1 J 

To this end, we use the values of the RG-invariant masses in the one-loop approximation, 
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FIG. 4: Left: The relative errors A.n[L] for N = 2, 3, and 4 of the truncated FAPT expansion, 
given by Eq. (44), with respect to the result of the resummation procedure, represented by Eq. 
(43). Right: The width T H _^. bb as a function of the Higgs-boson mass Mh in the resummed (blue 
solid line) and the truncated (at the order N) FAPT for three different values of N represented by 
dashed red lines. The short-dashed line corresponds to N = 1, the dashed one to N = 2, and the 
long-dashed line to N = 3. 

which we have extracted from two different sources: m b = 8.22 GeV [53] and rh b = 
8.53 GeV [54], with details being provided in Appendix E. 

In Fig. 4 we show the relative errors, given by (45), for N = 2, N = 3, and N = 4 in 
the probed range of L G [11.7, 13.6]. We see that already r^TjL; 2] gives correct results to 
better than 2.5%, whereas r^TJ-L; 3] reaches a still better accuracy at the level of 1%. This 
means that, on practical grounds, there is no need to calculate further corrections, because 
in order to be correct at the level of 1%, it is actually sufficient to take into account only the 
first three coefficients up to ds. This conclusion does not change if we vary the parameters 
of the model Pn{t)- To be more precise, varying the coefficients di in a reasonable range — in 
correspondence to their order, say, about 5% for d 2 up to 30% for d 5 — we induce changes 
of the parameter c on the level of about 5% which leave the main results (and conclusions) 
unchanged. By the same token, we can conclude that the quality of the convergence of 
the considered series in FAPT is quite high with a tolerance of only a few percent. Let us 

TABLE IV: Perturbation Theory (PT) coefficients d n for the Ds series in the deformed versions of 
the model in comparison with the original one. 



PT coefficients 


di d 2 dz di d 5 


pQCD results with N f = 5 [32] 


1 7.42 62.3 620 


Normal Model (40) with c = 2.43, 5 = -0.52 


1 7.50 61.1 625 7826 


Enhanced Deformation (40) with c = 2.62, 5 = —0.50 


1 7.85 68.5 752 10120 


Reduced Deformation (40) with c = 2.25, 5 = -0.51 


1 6.89 52.0 492 5707 



now expand our statements about the uncertainties of our results with regard to the model 
generating function P(t). To this end, we deform our original model with c = 2.43 and 
5 = —0.52 in order to enhance or reduce the magnitude of the last known coefficient d^ and 
the value of the still uncalculated coefficient d^ (for details see Table IV). The results of this 
variation are shown in the left panel of Fig. 5 in the form of a strip, the upper boundary of 
which is formed by the enhanced version of the P(t) model, whereas the lower boundary of 



17 



the strip corresponds to the "reduced" version of the model. We see that the uncertainty 
caused by this deformation is less than 0.5%. 




FIG. 5: Left: The width r^_^ as a function of the Higgs-boson mass Mjj in the resummed FAPT 
expansion, using different models for the generating function P(t) — see Table IV. Right: The 
width r^? 6 £ as a function of the Higgs-boson mass Mjj in the resummed FAPT, varying both 
the generating function P(t) in accordance with Table IV and also the mass ihrn with 5rhn-\ = 
±0.1 GeV. In both panels the upper strip corresponds to the Penin-Steinhauser estimate mm = 
8.53 GeV [54], whereas the lower one derives from the value mm = 8.21 GeV determined by Kiihn 
and Steinhauser in Ref. [53]. Here and also in Fig. 6 we indicated on the abscissa the window of 
the mass values of the Higgs-boson still accessible to experiment. 

Now we want to discuss what changes are induced when one applies for the same kind 
of analysis the FAPT resummation approach at the two-loop order. In that case, there is a 
technical complication: The evolution factors are not simple powers a"[L], but more involved 
expressions like a^°[L] [1 + C\ a]" 1 , as one sees by comparing with Eq. (D5). For this reason, 
the result of the resummation is more complicated, finally amounting to Eq. (D17). 

Here we resort only to graphical illustrations of our results. In the left panel of Fig. 6, we 
discuss the convergence properties of the decay widths, truncated at the order N, relative to 
the resummed two-loop result r^_^. From this, we infer that our conclusions drawn from 
the one-loop analysis remain valid. Indeed, T^^r[L; 2] is correct to better than 2%, whereas 
r^^fL; 3] reaches an even higher precision level of the order of 0.7%. 

In the right panel of Fig. 6, we show the results for the decay width V™ ^(Mjj) in the re- 
summed two-loop FAPT, in the window of the Higgs mass allowed by existing experiments — 
LEP and Tevatron [55]. Comparing this outcome with the one- loop result — upper strip in 
the same panel of this figure — reveals a 5% reduction of the two-loop estimate. This re- 
duction consists of two parts: one part (~ +7%) is due to the difference in the mass rh in 
both approximations, while the other (« —2%) comes from the difference in the values of 
Rs(Mh) in the one-loop and the two- loop approximations. 

In our predictions we have considered two different options for the values of the RG- 
effective 6-quark mass, mb(m%), which we have taken from two independent analyses. One 
value originates from Ref. [54] (mb(mf ) = 4.35 ± 0.07 GeV), while the other was derived 
in Ref. [53] yielding a somewhat deviating result, notably, m^m^) = 4.19 ± 0.05 GeV. 
In Fig. 6 we show all two-loop quantities adopting for the RG-effective 6-quark mass the 
result obtained by Penin and Steinhauser in Ref. [54]. The relative difference between these 
two choices is of the order of 4%, so that squared masses differ by 7%. This means that the 
corresponding curves for the Kiihn-Steinhauser value of m&(rnf ) [53] can be actually obtained 
from those shown by downsizing them with an overall reduction factor of 7%. From this we 
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FIG. 6: Left panel: The two-loop width T H ^ b i as a function of the Higgs boson mass Mjj in the 
resummed (black solid line) and the truncated (at the order N) FAPT. Here, the short-dashed 
line corresponds to N = 1, the dashed one to N = 2, and the long-dashed line to N = 3. Right 
panel: The two-loop width ^ as a function of the Higgs-boson mass Mu in the resummed 
FAPT is shown (lower strip), varying the mass in the interval rh{2) = 8.22 ± 0.13 GeV according 
to the Penin-Steinhauser estimate mf,(m|) = 4.35 ± 0.07 GeV [54]. The upper strip shows the 
corresponding one-loop result. The window of the mass values of the Higgs-boson which is still 
accessible to experiment is explicitly indicated. 



conclude that the real theoretical uncertainty of the Higgs-boson width in this decay channel 
is de facto determined by the upper boundary of the Penin-Steinhauser estimate [54], with 
the lower boundary of the Kuhn-Steinhauser [53] being in the range of 6% (6=(5+7)/2). 



V. ADLER FUNCTION OF THE VECTOR CORRELATOR AND R e+e -^ hadrons 

So far, we have discussed only the Adler function related to the scalar correlator. But 
this sort of considerations can be applied to the vector correlator as well. To be specific, 
we are interested in modeling the generating function of the perturbative coefficients d n (see 
the first row in Table I) of the Adler function of the vector correlator (labeled below by the 
symbol V) [33, 34] 



OO / 

Dy[L] = d + di ^2 I 

n=l ^ 



ot s [L] 



7T 



(46) 



To account for the n-dependence of the coefficients d n , in accordance with the asymptotic 
model of (34), we write 

d V n=c n - l6 -^^T(n) } (47a) 
which can be derived from the generating function 

P (t) - Se-^-(t/c)e-^ 

Mt) ~ c{P^Y) ■ (47b) 

Our predictions, obtained with this generating function by fitting the two known coefficients 
d 2 and d 3 and using the model (47), have been included in Table V. 5 We observe a good agree- 



5 Note that is automatically equal to unity. 
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ment between our estimate dj = 27.1 and the value 27.4, calculated recently by Chetyrkin et 
al. in Ref. [33, 34]. Would we use instead a fitting procedure, which would take into account 
the fourth-order coefficient d<± in order to predict d$, we would have to readjust the model pa- 
rameters in (47) to the new values {c = 3.5548, 5 = 1.32448} -> {c = 3.5526, 5 = 1.32453}. 
These findings provide evident support for our model evaluation, and we may expect that 
our procedure will work in other cases as well. 

In order to explore to what extent the exact knowledge of the higher-order coefficients 
d n is important, we employed our model (47) with different values of the parameters c and 
5: c = 3.63 and 5 = 1.3231. These values are, roughly speaking, tantamount to replacing 
the exact value of the coefficient d^ = 27. 4 by something approximately equal to the NNA 
prediction obtained in [48, 52]. The difference between the analytic sums of the two models 
in the region corresponding to Nf = 4 is indeed very small, reaching just a mere 0.2%. 



TABLE V: Coefficients d n for the Adler-function series with Nf = 4. The numbers in the square 
brackets denote the lower and the upper limits of the INNA estimates. 





PT coefficients 


di 


d 2 


d 3 


di 


d 5 


1 


pQCD results with Nf = 4 [33, 34] 


1 


1.52 


2.59 


27.4 






Models and estimates 




2 


Model (47) with c = 3.555, 6 = 1.3245 


1 


1.52 


2.59 


27.1 


2024 


3 


Model (47) with c = 3.553, 6 = 1.3245 


1 


1.52 


2.60 


27.3 


2025 


4 


Model (47) with c = 3.630, 8 = 1.3231 


1 


1.53 


2.26 


20.7 


2020 


5 


"NNA" prediction of [48, 52] 


1 


1.44 


13.47 


19.7 


579 


6 


"INNA" prediction of App. C 


1 


1.44 


[3.5,9.6] 


[20.4,48.1] 


[674, 2786] 


7 


"FAC" prediction of [36, 43] 








8.4 ±18 


152 



Then, we have 



N d v 



(48) 
(49) 



n=l 



while the global- APT resummation result for Dy PT [L] is given in Appendix B by Eq. (B2b). 
The relative errors 



K(Q 2 ) = i 



Z^ PT [ln(Q 2 /A 2 );iV] 
Z^ PT [ln(Q7A2)] 



(50) 



evaluated in the range Q 2 G [2, 20] GeV 2 for three values of the expansion = 1,2, 3, are 
displayed in Fig. 7. We observe from this figure that already _Dy PT [L; 1] provides an accuracy 
in the vicinity of 1%, whereas _Dy PT [L; 2] is smaller then 0.1% in the interval Q 2 = 1—20 GeV 2 . 
This means that there is no real need to calculate further corrections. Staying at the level of 
being correct to a better accuracy than 1%, it is virtually enough to take into account only 
the terms up to d 2 - This conclusion is quite robust against the variation of the parameters 
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FIG. 7: The relative errors A^(Q 2 ) evaluated for different values of N: N = 1 (short-dashed red 
line), N = 2 (solid blue line), and N = 3 (dashed blue line) of the truncated APT given by Eq. 
(49), in comparison with the exact result of the resummation procedure represented by Eq. (B2a). 

of the model Py(t). The main outcome here may look somewhat surprising: In fact, the 
best order of truncation of the FAPT series in the region Q 2 = 2 — 20 GeV 2 is reached by 
employing the N 2 LO approximation, i.e., by keeping just the c?2-term. 

VI. CONCLUSIONS 

In this work we have given considered in detail relations among popular perturbative 
approaches: FOPT, CIPT, and FAPT. We proved that in the Minkowski region both CIPT 
and FAPT produce for the i?-ratios which are related to the corresponding Adler functions — 
see Eqs. (22)-(24) — the same results. These results do not coincide — for any fixed order of 
the perturbative expansion — with those obtained with FOPT. 

We also considered in detail the resummation properties of non-power-series expansions 
within FAPT. In particular, we have given analysis of the Adler function of a scalar, D$, and 
a vector, Dy, correlator presenting results at the two- loop running of the coupling. Using 
a particular generating function for the coefficients of the perturbative expansion, which 
embodies information about their asymptotic behavior, we derived results for the whole 
series by resumming it. We used this key feature of the non-power-series expansion within 
FAPT in order to reduce the theoretical uncertainties in obtaining estimates for crucial 
observables, like the decay width of the Higgs boson into a bb pair, relevant for the Higgs 
search at the Tevatron and the LHC (Sec. IV). 

Employing an appropriate generating function, we estimated the values of the coefficients 
d% and and found that they are close to those computed by other groups using different 
perturbative methods. Moreover, we were able to resum the whole series and extract reli- 
able predictions for the width T H ^. b i in terms of the Higgs-boson mass in comparison with 
analogous estimates at fixed orders N = 1, 2, 3 of the perturbative expansion. This allowed 
us to gauge the accuracy of the truncation procedure of the FAPT expansion, estimating the 
relative errors for N = 2, 3, 4. We found that the convergence of the FAPT non-power-series 
expansion involves uncertainties at the level of only a few percent. A cautious conclusion 
from this is that the reached accuracy of the order of 1% at the truncation level of N 3 LO is 
comparable with, or even slightly better than, the 2% uncertainty involved in the m^ijn 2 ) 
estimates. For this reason, there is no real need to take into account higher-order corrections. 
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Our Higgs-decay predictions are given in graphical form in Figs. 4 (one-loop-order run- 
ning), 5 (one- loop-order running), and 6 (two-loop-order running). In Fig. 5 we pre- 
sented our one-loop resummed FAPT results for two different estimates of the mass scales: 
77i(i) = 8.53 GeV [54] (upper strip) and 777,(1) = 8.21 GeV [53] (lower strip). The graphi- 
cal representation of the two-loop decay width in the resummed FAPT is displayed in Fig. 
6 — lower strip — in comparison with the one- loop result (upper strip). In this calculation we 
varied the mass in the interval 777,(2) — 8.22 ±0.13 GeV, implementing the Penin-Steinhauser 
estimate m b (ml) = 4.35 ± 0.07 GeV [54]. The corresponding predictions for the other op- 
tion, offered by the Kuhn-Steinhauser estimate m b (ml) [53], can be obtained by reducing 
the previous results by an overall factor of about 7%. On the experimental side, one should 
keep in mind that the decay mode H — > bb is very challenging, but due to be measured by 
the ATLAS Collaboration at LHC [56]. 

In Sec. V we turned our attention to the Adler function, related to a vector correlator, in 
order to perform the FAPT resummation in the same manner as we did for the Higgs-boson- 
decay width. We used the generating function (47b) which depends on two parameters. 
To validate the robustness of our predictions, we varied the values of these parameters and 
found that this variation exerts only a small effect on the predicted value of d 5 and on the 
resummation result. It turns out that for a fixed number of flavors Nf = 4, the obtained 
results are within the limits set by the INNA method, being, however, incompatible with 
both the NNA prediction of [48, 52] and also the FAC one [36, 43]. We also estimated the 
influence of heavy-quark thresholds crossing in Appendix B. 

Bottom line: We provided evidence in terms of two concrete examples that FAPT can 
provide accurate and robust estimates for relevant observables that are otherwise inaccessi- 
ble by FOPT or CIPT. Moreover, using the suggested resummation approach within FAPT, 
we could optimize the truncation of the perturbative non-power-series expansion, thus min- 
imizing the truncation uncertainty. 
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Appendix A: Relations between R and D 



In this Appendix we derive expression (6) basing our reasoning on the useful relations 
between R (LHS) and D (RHS) 




(Alb) 



(Ala) 
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worked out in [50]. Equation (Alb) can be rewritten in the form of a dispersion relation 
between two specific quantities R and D to read 



D 



R 



Q 2 



ds 



-R = Q 



f 

Jo 



ds 



(s + Q 2 ¥ 



sin(7r5) f fi 
ty5 



,2\ * 



2 \ & 



From this equation we deduce the following correspondence between the terms R (LHS) and 
D (RHS): 



R = 



sm(n5) f /i 
ty5 



2\ S 



Q 2 



2\S 



D. 



(A3) 



The powers of ln(<5 2 / /x 2 ) in D(Q 2 , /i 2 ) generate, by means of R [D] in Eq. (4), the 7r 2 -terms 
in R(s,fi 2 ). To demonstrate how this happens, we take 2k times derivatives with respect to 
the variable n5 on both sides of Eq. (A3) and set 5 = 0. Then, we obtain a new relation 
between D and R, given by 



D (ln(g 2 //i 2 )) 



2k 



("I) 



2k + 1 



+ ... , 



(A4) 



where the ellipsis on the RHS indicates that the terms with the powers of ln(s//i 2 ) in R 
should, ultimately, disappear upon setting s = /i 2 , because only the even powers will survive. 
This can be traced back to the fact that the (ln(<5 2 //i 2 )) m terms in D originate from the 
expansion of the RHS of equation 

d n a n s (Q 2 ) =d n a n s ^ 2 ) (l + a s (/i 2 )/3 ln(g 2 //i 2 ))~ n 

that leads to the expression 



dn <(Q 2 ) = d n a n s J2 t !^ ! , Mo HQ 2 I^)T ■ 



m=0 



(n- 1)! (m)\ 



(A5) 



Finally, substituting relation (A4) into the expansion entering Eq. (A5), one arrives at Eq. 
(6) 



d n a n s ->■ 



(n-l + 2k)\ , k {a s M 2k 
^{^(n- 1)! (2A;)! ( ij 2k + 1 



(A6) 



Appendix B: Complete global expressions for D and R 

In this appendix, we derive explicit expressions for i? slob [L] and D glob [L] valid in the global 
scheme of APT. The spectral density p^ ob (<r) is related to the spectral density p n (cr; Nf) [7, 
57, 58]: 

6 

fC » = 9[a < ml] p n (a; 3) + ^ [mj < a < m 2 f+1 ] p n (a; N f ) , (Bl) 

where — in order to make our formulae more compact — we wrote m 4 = m c , m 5 = m&, m 6 = 
m t , and m 7 = +oo. 
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Then, as was shown in [27, 29], we have for 



D FAPT [L] 



R^[L] = d + d 1 ^e(L f <L<L f+1 )^i 1 [L + \ f -j--J 

5 6 

+ d 1 J2HL f <L<L f+1 ) 

f=3 k=f+l 

6 // rL 

f=3 \Y II r 




'+^1 [L a + X f ;N f ] dL a 



I _|_ e L-L a -t/P f 



+ d 1 ]T«A / [ J M]» p , 

/=4 



(B2a) 



(B2b) 



where A/ = ln(A|/Aj) describes the shift of the logarithmic argument due to the change of 
the QCD scale parameter A3 — > A/, (with Lf = ln(m 2 f / Af)) , L3 = —00, and L-j = +00. The 
second term in Eq. (B2a) presents a natural generalization of the fixed- Nf formula (30) for 
the case of using different QCD scales Aj for each fixed- Nf interval. The last term in Eq. 
(B2b) appears because of the continuation of R gloh [L] at the heavy-quark thresholds: 



Af^t}^^ Lf + Xf - —; f 

Pi ■ 



L f + A/_i 



fif- 



-J 



(B3) 



In the Euclidean domain, the corrections ((Af[L, t])) p to the naive expectation formula are 
defined by 



A f [L,t] 



1 p 1 [Lf + Xf-tx/Pf-Nf] t 



Pi [1 + 



yL—Lr—tx 
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dx 



-I 



1 p 1 [L f + X f - 1 - tx/pf-i, Nf.!] t 
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In contrast to the Minkowski case, they explicitly depend on L. 

Consider now the extension of these summation techniques to global FAPT, i.e., when 
one takes into account heavy-quark thresholds. To be more precise, we will deal with the 
summation of the following series: 



Dl 

DFAPT 
1 L u 



[L] 
[L] 



d 



Arm 

Vlt ob [L] 



+ d L dr. 



n=l 



•A-n+ v [L] 
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Note that due to the different relative normalization of ^4 n+i ,[L] (2l n+J ,[L]) and ^4^. b „[L] 
(^■n+ u [L]), the coefficients d n in Eqs. (27) and (B5) are also different. In order to obtain 
a generalization of the resummation procedure, given by Eq. (32), we propose to apply it 
to the spectral densities p^J^La]. This will be done for every A^-fixed integration region 
causing its reduction to p n+u [L a + A/; Nf]. Subsequently, application of Eq. (32) to sum up 
the series J2 p n +u[L a + X f ;Nf] ((r _1 )) p will yield 



n>l 



J2Pn + ALa + \ f ;Nf]((t n - 1 )) p = (U +u 



n>l 



L a + X f - 



Nf 



(B6) 
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As shown in [29], this procedure generates in the Minkowski region the following answer for 
the analytic image of the entire sum: 



< APT [L] = d Sir [L}+d 1 J^HL f <L<L f+1 )(^ 1+ ^L+X f ---J 

6 6 

+ d l J2HL f <L<L f+l )J2 «A fc a 1+ „M» A 

/=3 k=f+l 



(B7) 



with P v {t) taken from Eq. (33). 



Appendix C: Improved Naive Non-Abelization procedure 

Following the analysis of Ref. [28], we consider an expansion of the perturbative coeffi- 
cients d n in a power series in b , b±, . . ., i.e., 

d 3 = b 2 d 3 [2, 0] + h d 3 [0, 1] + b d 3 [l, 0] + d 3 [0, 0] , 

as opposed to the standard expansion in a power series in Nf (i.e., the number of flavors), 

d 3 = Nj d 3 (2) + Nj d 3 (l) + N° f d 3 (0). 

Here, the first argument no of the coefficients d n [no, ni, . . .] corresponds to the power of 
bo, whereas the second one n\ refers to the power of bi, etc. The coefficient d n [0, 0, . . . , 0] 
represents the "genuine" corrections which are associated with the coefficients fej in the power 
rii = 0. If all the arguments of the coefficient d n [. . . , m, 0, . . . , 0] to the right of the index m 
are equal to zero, then, for the sake of a simplified notation, we will omit these arguments and 
write instead d n [. . . , to]. Applying this terminology, we obtain the following representation 
for the next few coefficients 

d 3 = 6grf 3 [2]+6i d 3 [0,l] + b d 3 [l] + d 3 [0], (CI) 
d 4 = 6g d 4 [3] + h b d 4 [l, 1] + h d 4 [0, 0, 1] + b 2 d 4 [2] + h d 4 [0, 1] + b d 4 [l] + d 4 [0] , (C2) 
d 5 = bl d 5 [4] + 6i 6§ rf 5 [2, 1] + fcg dg [0, 2] + b b 2 d 5 [1 , 0, 1] + b 3 d 5 [0, 0, 0, 1] + O (bp . (C3) 

The same ordering of the /3-function elements applies with respect to the higher coefficients 
d n as well. Employing the standard NNA approximation, one estimates d n from the first term 
in the equations above, namely, d n ~ b^d n [n — 1]. We are going to improve this estimate by 

including other terms of the same order as 6q by means of the relation bi ~ 0(&q +1 ^). [For the 
reader's convenience, we have underlined the coefficients of these terms in Eqs. (C1)-(C3).] 
Note that in the MS scheme, the proportionality coefficients q = bi/b^ 1 ^ are 

Ci w 0.74 ; c 2 ~ 0.7 ; c 3 ps 1.82 ; c 4 = unknown . (C4) 

All these coefficients are of order unity and their values have been estimated for Nf = 
4, meaning that there is no reason to neglect them in Eqs. (C1)-(C3). This observation 
allows one to generalize the large bo- approximation, by taking into account all terms with 
the underlined coefficients in Eqs. (C1)-(C3) that belong to the same order in b . Within 
this setup, the element d n [n — 1] can be easily obtained from the term d n (n — 1), i.e., via 



25 



d n [n — 1] = (— |) n 1 d n (n — 1). Moreover, the coefficients d n (n — 1), related to the vector 
correlator, and the sum rules pertaining to deep-inelastic scattering were obtained for any n 
in [48, 52], respectively. 

The determination of the remaining underlined elements in Eqs. (C1)-(C3) is a difficult 
task that has been partially carried out in [28] for the particular case of the Adler D-function 
in the representation of Eq. (CI). Strictly speaking, it was found that d 3 [0, 1] — 0.4rf 3 [2] 
and d%[0, 1] ps 1.2, and it was proved that both elements are of order unity. With this in 
mind, we propose to use in the expansions of D$, [cf. (39)], and Dy, [cf. (46)], the following 
relation for the considered coefficients 

| rf ra [w ,7ii, . . .] [ = \d n [n - 1]| . (C5) 

Surprisingly, this rough approximation leads to reasonable results for the coefficients d n 
(especially when compared with those found with the NNA method), as one can see from 
the entries dubbed "INNA" in Tables I, III, V. As an illustrative example of this procedure, 
use Eq. (C5) in Eq. (C2) to predict the coefficients 

d 4 w t? 4 NNA = d 4 [3] + hb ± 62) + O(b 2 ) (C6) 
d 5 » 4 NNA = d 5 [4] (bt + hb 2 + 6? + 60 62 ± 63) + O(&o) , ( C7 ) 

with their values being given in Tables III and V. This kind of approximation is based 
upon the condition (C5) and is in line with the underlying assumptions of the original NNA 
procedure. 



Appendix D: Two-loop results 
1. Recurrence relations 

The expansion of the /3-function in the two-loop approximation is given by 

d (* S [L}\_ b J^V_ b faM\^ (D1) 



dL\ 4tt J " \ Air ) V 4vr 

where L = ln(/i 2 /A 2 ) and 

11 4 34 / 20 \ 

b = - C A - - T R N f , b ± = - C\ - [AC? + - C A j T R N f (D2) 

with C F = (A c 2 - 1) /2A C = 4/3, C A = N c = 3, T R = 1/2, and N f denoting the number 
of active flavors. Then, the corresponding two- loop equation for the coupling a = b a/(4ir) 
reads 

<h ' r2,[L] - " 2 [L](l + Cl a (2) [L]) wi,l,r ; (D3) 



dL - u (2)l _ JV . , -_ WL - Jy b . 

Still higher beta-function coefficients, e.g., b 2 , 63, can be found in [28, 59]. This equation 
immediately generates the following recurrence relation 

--_^-i=a^]+ciag 2 [L] (D4a) 
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for consecutive powers of the coupling constant. Due to its linearity, this relation remains 
valid also for the analytic images of the coupling's powers: 



1 d 
n dL 



(D4b) 



where, as it has already been used in Sec. Ill, J-[L\ denotes one of the analytic quantities 
«4.( 2 )[L], 2l( 2 )[L], or p^[L]. Quite analogously, we obtain the following generalization of this 
relation, pertaining to fractional coupling-constant indices, viz., 



-TT J~\ n+v\P\ — J 7 n+l+v[L] + Cl J 7 n +2+v[L] 



n + v dL 



(D4c) 



In the two-loop approximation we have a different evolution for the running mass, which 
reads 



= m 



1 + 



where 



4 7T 2 
(2) ^ fl (2)W ) 



4 7T 

[i + Cl a (2) (g 2 )] 



and rh{2) is the renormalization-group-invariant mass. 



(D5a) 
(D5b) 

(D6) 



2. Resummation in FAPT for fixed Nf 

Consider here the following power series with v > 0: 



n>\ 



(D7) 



noting that for v = we would obtain the corresponding two-loop APT expression. Multi- 
plying both sides of Eq. (D4c) with t n+u and summing over n from n = 1 to 00, we arrive 
at 

t 

-P J t' u W„[L; f; F] dt' = t"[(l + j) (m[L; t; F] - T 1+V [L]) - a T 2+V [L]] , 


where p = d/dL. Differentiating this equation with respect to t, we obtain 

j t {*" [(l + j) (W„[L; t; T\ - T l+V [L\) ~ d ^2+u[L]\ }+pt» W U [L; t;F} = (D8) 

with the initial condition W U [L; 0; J 7 ] = Fi +ll [L\. In order to solve this equation, we consider 
the following function 



X V [L; t] = (l + y) t^'^n+AL] ~ Cl t V F 2+v [L] , 



rt>2 



X v [L;0] = 0. 



(D9a) 
(D9b) 
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Our initial series W[L; t; J 7 ] is related to this function by the evident relation 



W V [L; t; F] = T 1+V [L] + X V [L; t] + T 2+U [L] . 



Ci+t 



Hence, we find 



^X v [L-t]+i>-^—X v [L-i\ = -pt»F l+u \L\-p C -±±—F 2+u \L\. 
at c\ + 1 ci + 1 



Finally, using the substitution 



X v [L;t]=e-* T ®X v [L;t], 



we get 



dX v [L;t] 
dt 



Ci+t 



That implies the relations 
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X u [L;t] = -t 1+ »Jz» 



dz 



dF 1+v [L + r(t z)] citz dT 2 + v \L + r(t z)\ 



dL 
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C\ + t z 
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dF 1+u [L + r(tz) -r(t)} | c x tz dT 2+v \L + r(t z) - r(t)} 



dL 



Ci+tz 



dL 



where we took into account that p = d/dL and where we used the initial condition (D9b). 
Making use of this solution in Eq. (Did), we obtain 



i 

W V [L; t; T] = T 1+V [L] + — ^— U F 2+U [L] - t [ z v dz 

ci + t I J 



dT 1+v [L + r(t z) - r{t)\ 



dL 



^ C\ + t z 

Note that this equation can be recast in the equivalent form 



^tz dT 2+l/ [L + r(tz) — r(t)] 



dL 



]} 



W u [L;t;T] = T 1+V [L]- 



t 2 f VA dF 1+u [L + r(tz)-T(t)} 



c 1 +t 



z v dz 



dL 



+ 



Ci t 
C!+t 



1 

{ T 2+V [L - r{t)\ V + uj z"' 1 T 2 + V \L + r(tz) - r{t)\ dz} , (D12) 



with <5o 5 „ being a Kronecker delta symbol. The reason for this recasting is that it is more 
appropriate for realizing the limit c\ — > and obtain this way the one-loop expressions, 
given by Eqs. (32)-(33). We observe that the analytic sum of the initial power series in t 
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can be represented by means of the analytic couplings J r i +I/ [L] and T 2 + V \L\. Setting v — 0, 
we obtain the corresponding two-loop APT expression 



Ci+tJ 



t W + ^-rWl^^^.^ (D13) 



dL 



C-L+t 



Hence, the sum expressed via equations like (32) can be recast in convolution form to 
read 6 



+ 



S v [L]=d F v [L]+d 1 ((T 1+V [L} + 



J^2+u[L — r(t)\ 5 0jl/ 

C\ ~\~ i // P 



(D14) 



— — Jzdz\ + —T 2+ AL + r(tz)-r(t)] 

o 



On the other hand, for the APT case with v — 0, expressions like (27) can be rewritten as 



S[L}=d + d 1 ((Fi[L]- 



t f dJ= x \L + r(f) - r(t)] Cl t 



C!+t 



1 



dL 



dt' + ^HL-rit)} )) . (D15) 



Quite analogously it is possible to resum also the following expression 
£v,u! [L] = d B u - Ul [L] + d\ — — - ^ p B n+ i +v . tV1 [L] , 



n=0 



(D16) 



in which B v - Vl [L] = Ae,m a ^ 2 )[L] (l + c i a (2)) UL [L] is the analytic image of the two-loop 
evolution factor (D5). Then, the FAPT result for the resummation of this series is 

Ci t 



S U . U1 [L] = d B v , Vl [L] + d l « B l+v . Vl [L] + 



■ B 2+ u-m l [L] 



d,t 2 



dz ■ 



z^ dB lm [L + r(tz)-r(t)\ 



(ci + t) 1 ^ J (ci + tz) 



dL 



di c\ t 



dz ■ 



t 2 z dB 2+V]Vl [L + r(t z) - r(t)\ 



(ci + t) 1 -^ J ~ (Ci + t z)* Ci + tZ 




dL 



+ 



di Ci t 2 



dz ■ 



Vi z 



{c 1 + t) 1 -^ J (d + tz) 1+u i 





B 2+u;Ul [L + r(t z) -r(t)\ 



(D17) 



Using the exact derivative expression 

dB 2+U]Ul [L + r{t z) - r{t)} 
dz 



t 2 z \ dB 2+v - v \L + r(t z) - r(t)\ 



Ci+t z 



dL 



Here we set S = D FAPT or R FAPT depending on the specific choice of T = A gloh or T = 2l glob . 
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and integrating by parts we can rewrite (D17) in a more convenient way to read 

l-z/i 



[L] = d B V . V1 [L] + di (( B X+V]V1 [L] + 5^ V+V1 1 



Cl 



Ci+t 



B 2+1/ , Vl [L-r{t)] 



(ci + f ) 



/ dz 



z^ dB 1+U]Vl [L + r{tz)-r{t)} 
(d + tz)^ dL 



di Ci t 



d 3 v + 1} , B 2+V , U1 [L + r{tz)- r{t)] 



{Ci+tz) 



(D18) 



3. Resummation in global FAPT 

The formalism developed in the text and the previous appendices can be applied to the 
global case using spectral densities, which correspond to a fixed number of active flavors Nf, 
and are defined in the integration intervals Lf < L a < Lf +1 . Indeed, employing T V \L\ = 
p u [L] and L = L a , we can first perform the resummation before carrying out the spectral 
integration over L a . 

Let us study these operations in some more detail within the Euclidean FAPT. In that 
case, the global spectral density Pn+ U { a ) nas the form specified in Eq. (Bl). It can be 
rewritten in an equivalent way to read 



/=3 



(D19) 



where we used L 3 = — oo and L 7 = +oo. Then, we have for the sum of the global power 
series (in t) the following expression 



n>l f=3 Jh f 



in which the following abbreviation was used: 



/+ 



[L a + \f,N f ;t] dL a 



1 + e 



L-L a 



(D20) 



p s r[L;N f ;t] = £ t n ~ *p n+v [L\ N f ] = ( -^Y + ' J] 



n>l / n >! 

We see that this series has the form of Eq. (D7) and is equal to 



t 

Jf 



n-l 



p n+v [L\ . (D21) 



PT m [L;N,;t] 



1 

Jf 



l+v 



T- — ■ 



(D22) 



where the last argument p^ of W u means that everywhere in Eq. (D12) one should substitute 
T[L] — > p^[L\. Hence, the initial sum 



(D23) 



n>l 
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can be resummed in the following final form 

V 



[L; 7V>] ee 



1 

ft 



f=?, JL f 
t 



(D24a) 
(D24b) 



Appendix E: Pole and RG-invariant masses of the bottom quark 

We start our considerations by writing down the relation [60, 61] between the pole mass 
of the b quark, m b , on one hand, and the value of the running mass at the scale //* = m b (fil), 
which we call m b (ml), on the other hand: 



m b (m 2 b ) 
m b 



8 cc, (ml) 
1 iA_J^ _ 13.234 

3 71 



a s (ml) 



-i 2 



7T 



(El) 



Note that some authors prefer another version of this relation in which the pole mass mj is 
used as argument in the numerator of the left-hand side: m b (ml) — see [60-63]. 7 

The QCD scales are fixed via the normalization of the strong coupling in the cor- 
responding one- or two-loop approximation at the Z-pole, employing the condition 

-R e + e -^hadrons(a s (^|)) = 1-039, suggested in [33], where 



Qt 

iZe+e-^hadronsftts) = 1 + — + 1.4097 

71 



(E2) 



TABLE VI: The values of the QCD scale parameters Af and 6-quark masses pertaining to [53]. 



Loop order 


A 3 [MeV] 


A 5 [MeV] 


m b (ml) [GeV] 


m b [GeV] 


m b [GeV] 


1-loop 


184 


115 


4.21 


4.69 


8.22 


2-loop 


346 


195 


4.19 


4.84 


7.89 



TABLE VII: The values of the 6-quark masses employing the estimates of the Penin-Steinhauser 
paper [54]. 



Loop order 


A 3 [MeV] 


A 5 [MeV] 


m b (ml) [GeV] 


m b [GeV] 


m b [GeV] 


1-loop 


184 


115 


4.35 


4.84 


8.53 


2-loop 


346 


195 


4.34 


5.00 


8.22 



Tables VI and VII show the results for m b obtained in this work. In parallel, we present 
results for the pole mass m b , using as input the estimates for m b (m b l ) derived by Kiihn and 
Steinhauser in Ref. [53] — Table VI. In Table VII we present analogous results, determining 
m b from the estimates for m b (ml) derived by Penin and Steinhauser [54]. For the sake of 



7 We wish to thank A. Kataev for attracting our attention to this point. 
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completeness, we also quote here the recent estimates m^ra^) = 4.163 ±0.016 GeV [64] and 
m b (ml) = 4.164 ± 0.025 GeV [65]. 
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